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SPACES OF FINITE ELEMENT DIFFERENTIAL FORMS 

DOUGLAS N. ARNOLD 

In memory of Enrico Magenes, in gratitude for his deep and elegant mathematics, 
which taught us, and his profound humanity, which inspired us. 

Abstract. We discuss the construction of finite element spaces of differential forms which satisfy 
the crucial assumptions of the finite element exterior calculus, namely that they can be assembled 
into subcomplexes of the de Rham complex which admit commuting projections. We present two 
families of spaces in the case of simplicial meshes, and two other families in the case of cubical 
meshes. We make use of the exterior calculus and the Koszul complex to define and understand the 
spaces. These tools allow us to treat a wide variety of situations, which are often treated separately, 
in a unified fashion. 



1. Introduction 

The gradient, curl, and divergence are the most fundamental operators of vector calculus, ap- 
pearing throughout the differential equations of mathematical physics and other applications. The 
finite element solution of such equations requires finite element subspaces of the natural Hilbert 
space domains of these operators, namely i7(curl), and i7(div). The construction of subspaces 
with desirable properties has been an active research topic for half a century. Exterior calculus 
provides a framework in which these fundamental operators and spaces are unified and generalized, 
and their properties and inter-relations clarified. Each of the operators is viewed as a particular 
case of the exterior derivative operator d — d!^ taking differential fc-forms on some domain C 
to differential [k + l)-forms. We regard as an unbounded operator between the Hilbert spaces 
L^A^ and L^A^+^ consisting of differential forms with coefficients. The domain of d^ is the 
Hilbert space 

(1.1) HK^ ^{u^ L^A^ I du G L^A^+i }, 

and all the d^ and their domains combine to form the de Rham complex 

^ HA^ A HA^ A • • • HA"" 0. 

Differential 0-forms and n-forms may be identified simply with functions on and differential 1- 
forms and (n — l)-forms may be identified with vector fields. In three dimensions, we may use these 
proxies to write the de Rham complex as 

^ grad ,x curl tt/ t \ div j-o 

^ - — > i7(curl) > i7(div) — > 0. 

The finite element exterior calculus (FEEC) is a theory developed in the last decade [H O [6] 
which enables the development and analysis of finite element spaces of differential forms. One 
major part of FEEC is carried out in the framework of Hilbert complexes, of which the de Rham 
complex is the most canonical example. One important outcome of FEEC is the realization that the 
finite dimensional subspaces A^ C HA^ used in Galerkin discretizations of a variety of differential 
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equations involving differential fc-forms should satisfy two basic assumptions, beyond the obvious 
requirement that the spaces have good approximation properties. The first assumption is that the 
subspaces form a subcomplex of the de Rham complex, i.e., that dA^ C A^^^. The second is that 
there exist projection operators tt^ from HA^ to which commute with d in the sense that the 
following diagram commutes: 



The second major part of FEEC, into which the present exposition falls, is concerned with the 
construction of specific finite element spaces A^ of differential forms. A special role is played by 
two families of finite element spaces V~A^{Th) and VrA^iTh)^ defined for any dimension n, any 
simplicial mesh 7/^, any polynomial degree r > 1, and any form degree < /c < n. Both these spaces 
are subspaces of Hh^{^). The V' spaces with increasing k and constant r form a subcomplex 
of de Rham complex which admits commuting projections. The same is true of the VrA^ family, 
except in that case the polynomial degree r decreases as the form degree k increases. 

We also discuss cubical meshes. In this case, there is a well-known family of elements, denoted 
by Q~ A^ in our notation, obtained by a tensor product construction. As for the V~A^ family, 
the Q-A^ spaces with constant degree r combine to form a de Rham subcomplex with commuting 
projections. We also discuss a recently discovered second family on cubical meshes, the 5^ A^ family 
of [3J. Like the VrA^ family, the de Rham subcomplexes for this family are obtained with decreasing 
degree. Moreover for large r, the dimSrA^{Th) is much smaller dimension than dimQ~A^. The 
finite element subspaces of H^, i7(curl), and i7(div) from this family in three dimensions are new. 

The remainder of the paper is organized as follows. In the next section we cover some preliminary 
material (which the more expert reader may wish to skip). We recall the construction of finite 
element spaces from spaces of shape functions and unisolvent degrees of freedom. To illustrate we 
discuss the Lagrange elements and carry out the proof of unisolvence in a manner that will guide 
our treatment of differential form spaces of higher degree. We also give a brief summary of those 
aspects of exterior calculus most relevant to us. In section [3] we discuss the two primary families 
of finite element spaces for differential forms on simplicial meshes mentioned above. A key role is 
played by the Koszul complex, which is introduced in this section. Then, in Theorem |3.5[ we give a 
proof of unisolvence for the V~ family which we believe to be simpler than has appeared heretofore 
(a similar proof could be given for the Vr family as well). In the final section we review the 
two families mentioned for cubical meshes, including a description, without proofs, of the recently 
discovered Sr family. 



2. Preliminaries 

2.1. The assembly of finite element spaces. Recalling the definition of a finite element space 
[TT] . we assume that the domain C is triangulated by finite elements, i.e., its closure is the 
union of a finite set Th of closed convex polyhedral elements with nonempty interiors such that the 
intersection of any two elements is either empty or is a common face of each of some dimension. We 
denote by A^(r) the set of faces of T of dimensions d, so, for example, Ao(T) is the set of vertices of 
T, and A^(T) is the singleton set whose only element is T. We also define A(r) = [jo<d<n ^d{T)^ 
the set of all faces of T. In this paper we consider the two cases of simplicial elements, in which 
each element T of the triangulation is an n-simplex, and cubical elements, in which element is an 
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n-box (i.e., the Cartesian product of n intervals). To define a finite element space C i7A^(il), 
we must supply, for each element T G 7^, 

(1) A finite dimensional space V{T)^ called the space of shape functions^ consisting of differential 
/c- forms on T with polynomial coefficients. The finite element space will consist of functions 
u which belong to the shape function spaces piecewise in the sense that u\t ^ ^(^) 

T ^ Th (allowing the possibility that u is multiply- valued on faces of dimension < n). 

(2) A set of functionals V{T) R, called the degrees of freedom^ which are unisolvent (i.e., 
which form a basis for the dual space 1/(T)*) and such that each degree of freedom is 
associated to a specific face of / G A(r). 

It is assumed that when two distinct elements Ti and T2 intersect in a common face /, the degrees 
of freedom of Ti and T2 which are associated to / are in a specific 1-to-l correspondence. If u 
is a function which belongs to the shape function spaces piecewise, then we say that the degrees 
of freedom are single-valued on u if whenever two elements Ti ^ T2 meet in a common face, 
then the corresponding degrees of freedom associated to the face take the same value on u\t-^ and 
u\t2^ respectively. With these ingredients, the finite element space A^ associated to the choice of 
triangulation 7/^, the shape function spaces V{T)^ and the degrees of freedom, is defined as the 
set of all fc-forms on Vt which belong to the shape function spaces piecewise and for which all the 
degrees of freedom are single- valued. 

The choice of the degrees of freedom associated to faces of dimension d < n determine the 
interelement continuity imposed on the finite element subspace. The use of degrees of freedom to 
specify the continuity, rather than imposing the continuity a priori in the definition of the finite 
element space, is of great practical significance in that it assures that the finite element space can 
be implemented efficiently. The dimension of the space is known (it is just the sum over the faces 
of the triangulation of the number of degrees of freedom associated to the face) and it depends only 
on the topology of the triangulation, not on the coordinates of the element vertices. Moreover, the 
degrees of freedom lead to a computable basis for A^ in which each basis element is associated to 
one degree of freedom. Further, the basis is local, in that the basis element for a degree of freedom 
associated to a face / is nonzero only on the elements that contain /. 

The finite element space so defined does not depend on the specific choice of degrees of freedom 
in y(T)*, but only on the span of the degrees of freedom associated to each face / of T, and we 
shall generally specify only the span, rather than a specific choice of basis for it. 

2.2. The Lagrange finite element family. To illustrate these definitions and motivate the 
constructions for differential forms, we consider the simplest example, the Lagrange family of finite 
element subspaces of = HhP. The Lagrange space, which we denote VrA.^{Th) in anticipation 
of its generalization below, is defined for any simplicial triangulation Th in and any polynomial 
degree r > 1. The shape function space is V{T) = Vr(T)^ the space of all polynomial functions 
on T of degree at most r. For a face / of T of dimension d^ the span of the associated degrees of 
freedom are the functionals 

(2.1) u E Vr{T) ^ j{tvfu)q, q G Vr-d-l{f). f G A(r). 

In interpreting this, we understand the space Vs{f) to be the space M of constants if / is 0- 
dimensional (a single vertex) and s >0. Also the space Vs{f) = if s < and / is arbitrary. The 
notation tifu denotes the trace of u on /, i.e., its restriction. Thus there is one degree of freedom 
associated to each vertex namely the evaluation functional u u{v). For r > 2 there are also 
degrees of freedom associated to the edges e of T, namely the moments of u on the edge of degree 
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at most r — 2: 



u\-^ J{tTeu)q, qeVr-2{e). 



' e 

For r > 3 there are degrees of freedom associated to the 2-faces, namely moments of degree at most 
r — 3, etc. This is often indicated in a degree of freedom diagram, hke that of Figure [TJ in which 
the number of symbols drawn in the interior of a face is equal to the number of degrees of freedom 
associated to the face. 




Figure 1. Degrees of freedom for the Lagrange quartic space V4,JSP in 3 dimensions. 

A requirement of the definition of a finite element space is that the degrees of freedom be 
unisolvent. We present the proof for Lagrange elements in detail, since it will guide us when it 
comes to verifying unisolvence for more complicated spaces. 

Theorem 2.1 (Unisolvence for the Lagrange elements). For any r > 1 and any n-simplex T , the 
degrees of freedom \2.1\ are unisolvent on V{T) — Vr{T). 



Proof. It suffices to verify, first, that the number of degrees of freedom proposed for T does not 
exceed dimy(T), and, second, that if all the degrees of freedom vanish when applied to some 



u G V{T)^ then = 0. For the first claim, we have by (2.1) that the total number of degrees of 
freedom is at most 

x:#A,(T)dimP._,_,(M^)=x:(;;!0(';,') = (";^') =dimP.(T), 

where the second equality is a binomial identity which comes from expanding in the equation 
(1 + xY^^{l + xy~^ = (1 + and comparing the coefficients of on both sides. 

We prove the second claim by induction on the dimension n, the case n = being trivial. Suppose 



that u G Vr{T) for some simplex T of dimension n and that all the degrees of freedom in (2.1) 
vanish. We wish to show that u vanishes. Let F G A^_i(T) be a facet of T, and consider iipu^ 
which is a polynomial function of at most degree r on the (n — l)-dimensional simplex F, i.e., it 
belongs to Vr{F). Moreover, if we replace T by F and uhj iiFu'vn ( |2.1| ), the resulting functionals 
vanish by assumption (using the obvious fact that iTfiiFU — iTfU for / C F C T). By induction 
we conclude that iipu vanishes on all the facets F of T. Therefore, u is divisible by the barycentric 
coordinate function which vanishes on F, and, since this holds for all facets, u — (fllLo '^^)^ 



some p G Vr-n-i{T). Taking f — T and q — pin (2.1) we conclude that 



L 



which implies that p vanishes on T, and so u does as well. □ 
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Let us note some features of the proof, which will be common to the unisolvence proofs for all 
of the finite element spaces we discuss here. After a dimension count to verify that the proposed 
degrees of freedom are correct in number, or at least no more than required, the proof proceeded 
by induction on the number of space dimensions. The inductive step relied on a trace property of 
the shape function space V{T) = Vr{T) for the family, namely that tTFV{T) C V{F). Moreover, 
it used a similar trace property for the degrees of freedom: if G V{Fy is a degree of freedom for 
then the puUback otip G ^(r)* is a degree of freedom for V{T). The induction reduced 
the unisolvence proof to verifying that if u E V{T)^ the space of functions in V{T) whose trace 
vanishes on the entire boundary, and if the interior degrees of freedom (those associated to T itself) 
of u vanish, then u itself vanishes, which we showed by explicit construction. 

Finally, we note that the continuity implied by the degrees of freedom is exactly what is required 
to insure that the Lagrange finite element space is contained in H^: 

(2.2) VrA-^{Th) = {u e H^{n)\u belongs to Vr{T) piecewise }. 

Indeed, a piecewise smooth function belongs to H^{^) if and only if its traces on faces are single- 
valued. Thus if a function in H^{ft) belongs piecewise to Vr{T), its traces are single- valued, so the 
degrees of freedom are single- valued, and the function belongs to VrJ^^iTh)- On the other hand, 
if the function belongs to VrA.^{Th), its traces on faces are single- valued, since, as we saw in the 
course of the unisolvence proof, they are determined by the degrees of freedom. Thus the function 
belongs to H^{Q). 

2.3. Exterior calculus. For the convenience of readers less familiar with differential forms and 
exterior calculus we now briefly review key definitions and properties. We begin with the space of 
algebraic fc-forms on V: Alt^ V = {L : 'R \ fc-linear, skew-symmetric }, where the multilinear 

form L is skew-symmetric, or alternating, if it changes sign under the interchange of any two of its 
arguments. The skew-symmetry condition is vacuous if /c < 2, so Alt^ F = and, by convention, 
Ait^ y = M. If uj is any fc-linear map M, then skwcj G Alt^ V where 

a 

with the sum taken over all the permutations of the integers 1 to k. The wedge product Alt^ V x 
Alt^ V Alt^+^ V is defined 

u Afi= I ^ \ skw(a; (g) /x), a; G Alt^ V, ji e Alt^ V. 

Let vi^ . . . ^Vji form a basis for V . Denoting by 

S(A;, n) = { (ai , . . . , (jfc) G I 1 < ai < • • • < (jfc < n }, 

an element of Alt^ V is completely determined by the values it assigns to the fc-tuples {vcr-^ , • • • , '^o-^)? 
(J G Sfc. Moreover, these values can be assigned arbitrarily. In fact, the fc-form /x^-i A • • • A ii^^^ 
where fii, . . . , fin is the dual basis to vi, . . . , Vn, takes the fc-tuple {vai^ • • • ? ^^ak) 1? other 
such /c-tuples to 0. Thus dim Alt^ V = (^), where n = dimV. 

We define differential forms on an arbitrary manifold, since we will be using them both when the 
manifold is a domain in M"^ and when it is the boundary of such a domain. A differential /c-form 
on a manifold is a map u which takes each point x G to an element Ux G Alt^T^fi, where 
is the tangent space to O at x. In other language, is a skew-symmetric covariant tensor 
field on ft of order k. In particular, a differential 0-form is just a real-valued function on ft and a 
differential 1-form is a covector field. In the case is a domain in R^, then each tangent space can 
be identified with R^, and a differential fc-form is simply a map ft Alt^R^. In this context, it is 
common to denote the dual basis to the canonical basis for R^ by dx^^ . . . , dx^^ so dx^ applied to a 
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vector V — {v^ ^ . . . , v'^) G is its fcth component . With this notation, an arbitrary differential 
fc-form can be written 

crGX;(/c,n) 

for some coefficients : ^ M. 

Three basic operations on differential forms are the exterior derivative, the form integral, and 
the pullback. The exterior derivative duo of a fc-form c<j is a (A: + l)-form. In the case of a domain 
in R^, it is given by the intuitive formula 

d(a^ rfx^i A • • • A rfx^M = V dx^ A dx""' A • • • A dx^'K 

It satisfies (in general) the identity d^^^od^ = and the Leibniz rule d{(jj/\ii) = (dcj) A/x+( — l)^Cc;A 
[dfi) if is a fc-form. 

The definition of the form integral requires that the manifold be oriented. In this case we can 
define G M for c<j an n-form with n — dimfi. The integral changes sign if the orientation of 
the manifold is reversed. 

Finally, if F : ^ fiMs a differentiable map, then the pullback F* takes a fc-form on Q.' to one 
on by 

(F*cj)^(i;i, ...^vk)^ UF^^^{dFxVi, . . ..dF^Vk), x G 1^, vi,...,Vk e T^l^. 
The pullback respects the operations of wedge product, exterior derivative, and form integral: 

F*(cc;A/x) = (F*cc;) A(FV), F\du) ^ d{F'' u) , [ F'^uj = [ uj, 

for (jj and differential forms on The last relation requires that F be a diffeomorphism of 
with Q.' which preserves orientation. 

An important special case of pullback is when F is the inclusion of a submanifold into a larger 
manifold In this case the pullback is the trace operator taking a fc-form on Q.' to a fc-form on 
the submanifold Vt. All these operations combine elegantly into Stokes' theorem, which says that, 
under minimal hypothesis on the smoothness of the differential (n — l)-form uj and the n-manifold 



I iiuj — I duo. 
JdQ Jn 



If V is an inner product space, then there is a natural inner product on Alt^ V. Thus for a 
Riemannian manifold, such as any manifold embedded in R^, the inner product {cJxif^x) G M is 
defined for any fc- forms cj, jn and any x ^ ft. An oriented Riemannian manifold also has a unique 
volume form, vol, a differential n-form which at each point assigns the value 1 to a positively 
oriented orthonormal basis for the tangent space at that point. For a subdomain of the volume 
form is the constant n-form with the value dx^ A • • • A dx^ at each point. Combining these notions, 
we see that on any oriented Riemannian manifold we may define the L^-inner product of fc-forms: 



J r2 



The space L^A^ is of course the space of fc-forms for which ||cL;||^2yvfe •= (cj, oo)j^2f^k < oc, and then 
HA^ is defined as in (1.1). 
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3. Families of finite element differential forms on simplicial meshes 

Our goal now is to create finite element subspaces of the spaces Hh^ which fit together to yield 
a subcomplex with commuting projections. In this section the spaces will be constructed for a 
simplicial a triangulation Th of the domain C M^. Thus, for a simplex T, we must specify a space 
V{T) of polynomial differential forms and a set of degrees of freedom for it. 

3.1. The polynomial space T^^A^. An obvious choice for V(T) is the space 

of a differential fc-forms with polynomial coefficients of degree at most r. It is easy to compute its 
dimension: 

(3.1) dimP.A'^(r) = #E(fc,n)xdimP.(r)= + - ([^ ! ;) ^ ') ■ 

Note that dVr^^ C P^-iA^^^, i.e., the exterior derivative lowers the polynomial degree at the 
same time as it raises the form degree. Therefore, for each r we have a subcomplex of the de Rham 
complex: 

(3.2) Vr^^^ A n-lA^ A • • • A Pr-nA^ ^ 0. 

This complex is exact (we have left off the initial since the first map, d — grad acting on P^A^ 
has a 1-dimensional kernel, consisting of the constant functions). That is, if G PsA^ and duo — {) 



then uj — dfjL for some ji G T^^+iA^"^. We prove this in Corollary 3.2 below, using an elementary 
but powerful tool called the Koszul complex. The same tool will also be used to define the degrees 
of freedom for Vr^^iT)^ and to define an alternative space of shape functions. 

3.2. The Koszul complex. For a domain in C (but not a general manifold), the identity 
map may be viewed as a vector field. It assigns to an arbitrary point x G C the point itself 
viewed as a vector in and so an element of the tangent space T^fi. Contracting a fc-form uj with 
this identity vector field gives a (fc — l)-form KjUj: 

{k.u)x{vi, . . . , Vk-i) = Uj,{x, vi,..., Vk-i), X E Jl, vi,..., Vk-1 E 

Since ojx is skew-symmetric, kkoj = 0, that is, is a differential. It satisfies a Leibniz rule: 

k{cj a /x) = (kcj) a /X + (-1)^^ A (/^/i), 

for a fc-form cj and a second form fi. In particular hi{fuj) — fi^uj if / is a function. Also /^rfx^ = x\ 
These properties fully determine Thus 

n^dx"" A dx^) = dx^ - x^ dx\ K{dx'^ A dx^ A dx^) = dx^ A dx^ - x^ dx'^ A dx^ + x^ dx'^ A dx\ 

and so forth. If we identify 1-forms with vector fields, then corresponds to the dot product of the 
vector field with x (or, more properly, with the identity vector field). On 2-forms in 3-D, is the 
cross product with x, and on 3-forms it is the product of a scalar field with x to get a vector field. 

The Koszul differential k maps the space VrA.^ of differential fc-forms with coefficients in Vr{^) 
to P^+iA^~^, exactly the reverse of d. Thus both nd and dn map P^A^ to itself. The following 
theorem points to an intimate relation between and d, called the homotopy formula. In it we 
write T-LrA.^ for the fc-forms with homogeneous polynomial coefficients of degree r. 

Theorem 3.1 (Homotopy formula). 

{nd + dhi)oj = (A; + r)u, uj E 
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Remarks on the proof. The case A: = is Euler's identity x • gradj) = rp for p a, homogeneous 
polynomial of degree r. Using it, we can verify the theorem by direct computation. Alternatively, 
one may use Cartan's homotopy formula from differential geometry. For details on both proofs, see 
Theorem 3.1 of [5j. □ 



Corollary 3.2. The polynomial de Rham complex (|3.2|) and the Koszul complex 
are both exact. 



Proof. For the de Rham complex, it suffices to establish exactness of the homogeneous polynomial 
de Rham complex 

since then we can then just sum to get the result. We must show that if G T^sA^ and duo — 
then uj is in the range of d. Indeed, by the homotopy formula 

= (5 + k)~^{dHi + Hid)(jj = (5 + k)~^dHi(jj. 

A similar proof holds for the Koszul complex. □ 

Another important consequence is a direct sum decomposition: 

Corollary 3.3. For r > 1, < k < n, 

(3.3) HrA^ = ^Hr-iA^^^ e dnr+iA^~\ 

Proof. By the homotopy formula, any element of T-LrA^ belongs to kT-Lt-iA^^^ + dT-Lr^iA^~^ . More- 
over the intersection of these two spaces is zero, since if u belongs to the intersection, then du = 0, 
A^o; = 0, so = by the homotopy formula. □ 

3.3. The polynomial space V~A^. We now define a second space of polynomial differential forms 
which can be used as shape functions. We have 

VrA^ = Vr-lA^ e UrA^ = Vr-lA^ ^Ur-lA^^^ dUr+lA^'^ 

If we drop the last summand, we get a space intermediate between Vr-iA^ and VrA^: 

(3.4) V-A^ := Vr-iA^ + t^Ur-iA^^^ . 

Note that V~ A^ = VrA^ and P^A^ = Vr-iA"", but for < A: < n, V~ A^ is contained strictly 
between Vr-iA^ and VrA^. We may compute the dimension of kT-LtA^^ using the exactness of the 
Koszul complex and induction (see ^ Theorem 3.3]). This then yields a formula for the dimension 
of P-A^: 

T ^-xk fn^r\fr-\-k-V 
^^"^^^^^ n-J k 



Comparing this with (3.1), we have 

dim7^-A^ = dimP^A^ 
r + A: 

(showing again that the spaces coincide for 0- forms). 
Now 

dV'A^ C dVrA^ C Vr-iA^^^ c V-A^^\ 
so we obtain another subcomplex of the de Rham complex: 

(3.5) v-A^ A v;a^ p-A^ ^ 0. 
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Note that, in contrast to (3.2), in this complex the degree r is held constant. However, like (3.2), 
the complex (3.5) is exact. Indeed, 

where the penultimate equality follows from Corollary |3. 2| and the last equality is a consequence of 



the definition (|3.4|) and the homotopy formula Theorem 3.1 



3.4. The V~ A^(Th) family of finite element diff*erential forms. Let r>l,0<A:<n, and 
let 7/i be a simplicial mesh of C M^. We define a finite element subspace V~ A^iTh) of HA^ip). 
As shape functions on a simplex T e % we take V{T) = ■p-A''{T). As de errees of freedom we take 



(3.6) 



u G V-A^{T) ^ / {tif u)Aq, qe P.+^-^-iA^-^(/), / G A^(r), d>k. 



Note that, in the case A: = 0, ^(r) = Vr{T) and ([3^ coincides with ([2l]), so the space V~A^{Th) 
generalizes the Lagrange finite elements to differential forms of arbitrary form degree. We shall 
prove unisolvence for arbitrary polynomial degree, form degree, and space dimension at once. The 
proof will use the following lemma, which is proved via a simple construction using barycentric 
coordinates. 

Lemma 3.4. Let r > 1, < k < n, and let T he an n-simplex. If u ^ Vr-iA^{T) and 
(3.7) j^uAq^O, qeVr+k-n-i^'^-HT), 

then u = 0. 

Proof Any element of Vr-iA^(T) can be written in terms of barycentric coordinates as 

^ UcrdXcr^ A- • • AdXcr^, ^ Vr-l{T) . 

Now let 1 < z < n, and consider the trace of u on the face given by = 0. By the assumption 
that u G Vr-iA^{T), the trace vanishes. This implies that divides for any a G Ti{k^n) whose 
range does not contain i. Thus 

Ua = PaKl • • • for some p^ G Vr+k-n-liT), 

where a* G — k^n) is the increasing sequence complementary to a. Thus 

^ p^A^* • • • A^*_^dA^-L A • • • A dA^^, Pa ^Vr^k-n-i{T)- 

a^Ti{k,n) 



Choosing 



in (3.7), we get 







aeT.{k,n) 

= J uAq = J ^ p^A^* • • • A^*_^dAi A • • • A dA^. 



aeT.{k,n) 

However, the A^ are positive on the interior of T and the n-form dXi A • • • AdXn is a nonzero multiple 
of the volume form. Thus each pcr must vanish, and so u vanishes. □ 

Theorem 3.5 (Unisolvence for V~A^{Th))> For any r > 1, < k < and n-simplex T , the 
degrees of freedom (|3.6[) are unisolvent for V{T) — V~A^{T). 
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Proof. First we do the dimension count. The number of degrees of freedom is at most 



d>k d>k 



n-\-l\fr-\-k — l\fd 



j>0 



n+1 \ /r + A: - 1\ /j + A: 
j + k + l)\ j + k )\ j 



Simphfying with the binomial identities, 

fb\ fa\ fa — c\ ^-^ ( a \ ( c\ fa + c 



bj\cj \c)\a — h)^ \b+jj\jj \a — b 



the right-hand side becomes 



r + n\/r + /c-lx ^-,k^m\ 



It remains to show that if u ^ A^(T) and the degrees of freedom in ( |3.6[ ) vanish, then u 
vanishes. Since tTfV~A^{T) = P~A^(/), we may use induction on dimension to conclude that 
tr fU vanishes on each facet so u G V-A^{T). Therefore du G Vr-iA^^^{T). Moreover, 

/ duAp = ± / uAdp^O, p eVr+k-n^''~^~^{T), 
Jt Jt 

where the first equality comes from Stoke's theorem and the Leibniz rule, and the second from 
the hypothesis that the degrees of freedom for u vanish. We may now apply the lemma (with k 
replaced by /c + 1) to du to conclude that du vanishes. But the homotopy formula implies that for 
u G V~A^ with du — u ^ Vr-iA^. Using the interior degrees of freedom from (3.6), we may 
apply the lemma to u^ to conclude that u vanishes. □ 

It is easy to check that the degrees of freedom imply single- valuedness of the traces of elements 
of V~ A^(Th)^ so that they indeed belong to HA^. Moreover, it is easy to see that the complex (3.5) 
involving the shape functions, leads to a finite element subcomplex of the L? de Rham complex on 

V-A\Th) A V-A' i> (r^) . . . i> V-A^{Th). 

Using the degrees of freedom to define projection operators tt^ into V~ A^(Th) (the domain of tt^ 
consists of all continuous fc-forms in HA^{Q))^ we obtain projections that commute with d (this 
can be verified using Stokes' theorem), which is crucial to the analysis of the element via FEEC. 

3.5. The VrA^iTh) family of finite element differential forms. We may also use the full 
polynomial space VrA^(T) as shape functions for a finite element space. The corresponding degrees 
of freedom are 

(3.8) u G VrA\T) ^ [ {tTf u)Aq, g G / ^ MT). d>k. 

Note that in this case the degrees of freedom involve V~ spaces, defined through the Koszul complex. 
The analysis of these spaces is very parallel to that of the last subsection, and we will not carry it 
out here. Again, we obtain unisolvence, and a finite element subcomplex of the de Rham complex 

VrI^\Th) A Vr-lA\rh) A ... A Vr-nA^'iTh), 

which admits a commuting projection defined via the degrees of freedom. 
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3.6. Historical notes. In the case = 0, the two shape function spaces and P^A^ coincide, 

as do the spaces Vr-d-i^'^'^if) and V~_^A^d - k(f), f G Arf(T), entering ([3^ and ([3^. Thus 
the two finite element famihes coincide for 0-forms, and provide two distinct generahzations of the 
Lagrange elements to differential forms of higher degree. 

In n dimensions, n- forms may be viewed as scalar functions and the space HA^(Q) just cor- 
responds to L^(fi). The finite element subspace VrA^{Th) is simply the space of all piecewise 
polynomial functions of degree r, with no interelement continuity required. The space V~A^{Th) 
coincides with Vr-iA^{Th)- 

In two dimensions, the remaining spaces V~ h}{Th) and VrA.^{Th) can be identified, via vector 
proxies, with the Raviart-Thomas spaces [20j and the Brezzi-Douglas-Marini spaces [TOJ. In three 
dimensions, the V~ h}{Th) and V~ h?{Th) spaces are the finite element subspaces of i7(curl, Vt) and 
i7(div, fi), respectively, called the Nedelec edge and face elements of the first kind [18]. The spaces 
VrA.^{Th) and VrA.'^{Th) are the Nedelec edge and face elements of the second kind [T9]. Diagrams 
for the two-dimensional and three-dimensional elements are shown in Figures [2] and Figures [3j 

The lowest order spaces h^{Th) are very geometric, possessing precisely one degree of freedom 
per face of dimension fc, and no others (see the top rows of Figures [2] and |3j In fact these spaces 
first appeared in the geometry literature in the work of Whitney in 1957 [24j long before their first 
appearance as finite elements. In the 1970s, they were used by Dodziuk [13j and Dodziuk and Patodi 
|14| as a theoretical tool to approximate the eigenvalues of the Hodge Laplacian on a Riemannian 
manifold. This then played an essential role in Miiller's proof of the Ray-Singer conjecture |17| . 
The spaces VrA^{Th) also appeared in the geometry literature, introduced by Sullivan [2?^ "23]. In 
an early, largely overlooked paper bringing finite element analysis techniques to bear on geometry 
Baker [7] named these Sullivan- Whitney forms, and analyzed their convergence for the eigenvalue 
problem for the Hodge Laplacian. In 1988 Bossavit made the connection between Whitney forms 
and the mixed finite elements in use in electromagnetics [8j , in part inspired by the thesis of Kotiuga 
[16] . The first unified treatment of the V~A^ spaces, which was based on exterior calculus and 
included a unisolvence proof, was in a seminal paper of Hiptmair [15] in 1999. In the 2006 paper 
of Arnold, Falk, and Winther [5j, in which the term finite element exterior calculus first appeared, 
the Koszul complex was first applied to finite elements, simplifying many aspects and resulting in 
a simultaneous treatment of both the V~A^ and VrA^ spaces. 



4. Families of finite element differential forms on cubical meshes 

We now describe two families of spaces of finite element differential forms, which we denote 
Q~A^{Th) and SrA^{Th)^ defined for cubical meshes i.e., meshes in which each element is the 
Cartesian product of intervals. In some sense, the Q~A^ family can be seen as an analogue of the 
V~A^ family for simplicial meshes, and the ^^A^ family an analogue of the VrA^ family. The Q~A^ 
family can be constructed from the one-dimensional case by a tensor product construction, and is 
long known. By contrast, the SrA^ family first appeared in recent work of Arnold and Awanou 
[3j. Even in two and three dimensions, the spaces in this family were for the most part not known 
previously. 

4.1. The Q~A^ family. We describe this family only very briefly. A more detailed description will 
be included in a forthcoming study of the approximation properties of these spaces under non-affine 
mappings [4J . Suppose we are given a subcomplex of the de Rham complex on an element S C 
and a second such subcomplex on an element T C M^: 

V\S) A V\S) A • • • A 1/^(5), V\T) A V\T) A • • • A i/^(r). 
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fc = 1 



k^2 



r = 1 



r = 2 



r = 3 










fc = 1 







Legend: u - scalar (1 DOF); ^ - vector field (2 DOFs); 
- tangential component along edge (1 DOF). 



Figure 2. The V~A''{Th) and VrA^{Th) spaces in two dimensions. 
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VrA'': k = k = l k = 2 k = 3 




Legend: u - scalar (1 DOF); - tangential vector field (2 DOFs); - vector field (2 DOFs); 
- tangential component along edge or normal component on face (1 DOF). 

Figure 3. The V-A^{Th) and Vrh^{Th) spaces in three dimensions. 
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We may then construct a subcomplex of the de Rham complex on 5 x T by a tensor product 
construction which is known in the theory of differential forms; see, e.g., ^2T| p. 61]. The canonical 
projection ns S x T ^ S determines a pullback of z-forms on S to z-forms on 5 x T, so 7t'^V^{S) 
is a space of z-forms on S x T and, similarly, 7t^V^{T) is a space of j-forms on S x T. Thus we 
may define a space of fc-forms on 5 x T by 

i-\-j=k 

We take the space V^{S x T) as the shape functions for A:- forms on 5 x T. The construction of 
degrees of freedom for V^{S x T) is simple. If 77 G V^{SY is a degree of freedom associated to a 
face / of 5, and p E W(r)* is associated to a face g of T, we define 

7/ A p G [itIV\S) a ^^l/^'(r)]* C V^{S X T)*, 

by 

(7/ A p)(7r> A Ti^v) = r]{u)p{v), 

and associate the degree of freedom r] A p to f x g, which is a face of S x T. 

The Q~ family is defined by applying this tensor product repeatedly, starting with a finite element 
de Rham complex on an interval in one dimension. In one dimension the V~ and Vr de Rham 
subcomplexes coincide. On an interval /, the shape functions for 0- forms are V^{I) = Vr{I) 
with degrees of freedom at each end point, and moments of degree at most r — 1 in the interior. 
The shape function for 1-forms are V^{I) — Vr-i{I) with all degrees of freedom in the interior. 
Repeatedly using the tensor product construction just outlined, we obtain polynomial spaces and 
degrees of freedom on a box /i x • • • x C M^. We denote the shape function space so obtained 
by Q-K^{h X ••• y. In)- In n = 2 dimensions, for example, 

Q-A°(Ii X h) = Qrih X h) = Vr{h)®Vr{h), 
Q-K\h X h) = [Pr-lih) ® Vrih)] X [Vrih) ® Vr-lih)], A^h X h) = X h). 

Diagrams for these elements in two and three dimensions are shown in Figures |4j The space 
Q~A^{Th) is the standard Qr finite element subspace of H^{Q) and the space Q~A^{Th) is the 
discontinuous Qr-i subspace of L^(fi). The space Q~A^{Th) goes back to Raviart and Thomas [20j 
in two dimensions, and the Q~A^{Th) and Q~A'^{Th) were given by Nedelec in [18j. The spaces 
with r held fixed combine to create a finite element de Rham subcomplex, 

Q^A^Th) A Q-A\Th) A ... A Q-A^{Th), 

and the degrees of freedom determine commuting projections. 

Recently, Cockburn and Qiu [12j have published a different family of finite element spaces in two 
and three dimensions, that seems to be related to these. They begin with the complex formed by 
the full spaces QrA^, which lie between Q~A^ and Q~^iA^. That complex (which was discussed 
in [19]) does not admit commuting projections. Cockburn and Qiu define a small space of bubble 
functions that can be added to each of the spaces so that the resulting spaces remain inside Q~^iA^ 
but also form a de Rham subcomplex (with constant r) which admits commuting projections. 

4.2. A second family of finite element differential forms on cubes. The ^^A^ family pre- 
sented in this section was derived recently in [3J. It seems to be complementary to the Q~A^ family 
much as the VrA^ family complements the V~A^ family. To describe the new family we require 
some notation. A fc-form monomial in n variables is the product of an ordinary monomial and a 
simple alternator: 

m = (x^)^i • • • (x^)^- rfx^i A • • • A dx^'K 
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Figure 4. The Q~A^{Th) spaces in two and three dimensions. 
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where a is a multi-index and a G Tj{k^n). We define the degree of m to be the polynomial degree 
of its coefficient: degm = ai. The linear degree of m is more complicated: 

Idegm = #{z I = 1, ai ^ {ai,. . . ,ak}}, 

that is, the number of variables that enter the coefficient linearly, not counting the variables that 
enter the alternator. For example, if m = x^x^{x^)^ dx^^ then degm = 7, Idegm = 1. 

We now define the space of shape functions we shall use for fc-forms on an n-dimensional box, 
T. Viewing monomial forms as differential forms on T, we define Hri^^{T) C Hr^^{T) to be the 
span of all monomial fc-forms m such that degm = r and Idegm > /. Using this definition and the 
Koszul differential, we then define 

l>l 

Finally, we define the shape functions on T by 

defined for all r > 1, < A; < n. 

As the definition of the shape functions takes a while to absorb, we describe the spaces in more 
elementary terms in the case of three dimensions. 

• The space ^^A^, the polynomial shape functions for the space, consists of all polyno- 
mials u with superlinear degree sdegu < r. The superlinear degree of a monomial is its 
degree ignoring any variable that enters to the first power, and the superlinear degree of a 
polynomial is the maximum over its monomials. The criterion sdegi/ < r was introduced 
in [2] to generalize the serendipity elements from 2 to n-dimensions. 

• The space ^^A^, the shape functions for the i7(curl) space, consists of vector fields of the 
form 

with polynomials v'^, w'^, and u for which degv'^ < r, degw'^ < r — 1, sdegi/ < r + 1, and w'^ 
is independent of the variable x\ 

• The i7(div) space uses shape functions ^^A^, which are of the form 

v'^, v^) + curl(x^x^(it;^ - w^),x^x^{w^ - w'^), - w'^)), 

with degv'^ < r, degw'^ < r, and w'^ independent of the variable x\ 

• Finally the space ^^A^ simply coincides with Vr- 

In 21 we establish the following properties of these spaces (in any dimension): 

• degree property: 'P^A^(/^) C 5^A^(/^) C 'P^+^_/eA^(/^), 

• inclusion property: SrA^{I^) C 5r+iA^(/^), 

• trace property: For each face / of tif SrA^{I^) C 5^A^(/), 

• subcomplex property: d5^A^(/^) C 5^_iA^+i(/^). 

The degrees of freedom for SrA^{T) are quite simple: 

(4.1) u G SrA\T) ^ / {tTf u)Aq, qe Vr-2{d-k)A^~\f). f ^ Arf(T), d>k. 

These are illustrated in Figure [5| Notice that weighting function q is sought in a Vs space, not a 
Qs space. Moreover, as the face dimension d increases by 1, the degree s of the space used for q 
decreases by 2. A major result of ^ is a proof that the degrees of freedom are unisolvent. Further, 
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we show there that the resulting finite element spaces combine into de Rham subcomplexes with 
commuting projections: 

5,A0(r^) A Sr-lA\Th) A ■■■A Sr-nA^^iTh), 

in which the degrees r decrease, as for the Vr^^ iTh) spaces on simplices. 

For n- forms, the space Sr^^^{Th) is simply the discontinuous Vr space (but defined on boxes, 
rather than simplices). In 2-dimensions, the 0-form space SrJ^^{Th) is the well-known serendipity 
space, and the 1-form space is the rectangular BDM space defined in [lOj. Hence these spaces were 
all known in 2 dimensions. However, in 3 and more dimensions they were not. The 0-form space is 
the appropriate generalization of the serendipity space to higher dimensions, a space first defined 
in 2011 [2j. The space Srh? in 3-D is, we believe, the correct analogue of the BDM elements to 
cubical meshes. It has the same degrees of freedom as the space in [9J but the shape functions have 
better symmetry properties. For 1-forms in 3-D, 5rA^ is a finite element discretization of i7(curl). 
To the best of our knowledge, neither the degrees of freedom nor the shape functions for this space 
had been proposed previously. Finally, we note that the dimension of S~ A^{T) tends to be much 
smaller than that of Q~A^(r), especially for r large, as can be observed in Table [l] 



r r 
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49 
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12 


24 
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14 


22 


32 


44 


58 


2 


1 


4 


9 


16 


25 


36 




3 


6 


10 


15 


21 


28 
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8 


27 


64 


125 


216 


343 




8 


20 


32 


50 


74 


105 


1 


12 


54 


144 


300 


540 


882 




24 


48 


84 


135 


204 


294 


2 


6 


36 


108 


240 


450 


756 




18 


39 


72 


120 


186 


273 


3 


1 


8 


27 


64 


125 


216 




4 


10 


20 


35 


56 


84 














n 


= 4 

















16 


81 


256 


625 


1296 


2401 




16 


48 


80 


136 


216 


328 


1 


32 


216 


768 


2000 


4320 


8232 




64 


144 


272 


472 


768 


1188 


2 


24 


216 


864 


2400 


5400 


10584 




72 


168 


336 


606 


1014 


1602 


3 


8 


96 


432 


1280 


3000 


6048 




32 


84 


180 


340 


588 


952 


4 


1 


16 


81 


256 


625 


1296 




5 


15 


35 


70 


126 


210 



Table 1. Dimension of Q7A^(/^) and 5^A^(/^). 
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